Velocity Statistics Distinguish Quantum Turbulence from Classical Turbulence 
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By analyzing trajectories of solid hydrogen tracers, we find that the distributions of velocity in 
decaying quantum turbulence in superfluid 4 He are strongly non-Gaussian with l/v 3 power-law tails. 
These features differ from the near-Gaussian statistics of homogenous and isotropic turbulence of 
classical fluids. We examine the dynamics of many events of reconnection between quantized vortices 
and show by simple scaling arguments that they produce the observed power-law tails. 
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The pioneering work of Kolmogorov [JJ [S] remains the 
cornerstone of the statistical theory of classical incom- 
pressible turbulence. Kolmogorov made two key assump- 
tions: (1) local isotropy and homogeneity prevail, and 
(2) there exists an inertial range in which turbulent en- 
ergy is transferred from large to small scales indepen- 
dent of viscosity and generation mechanisms. Dimen- 
sional arguments then yield the spectral density E(k) = 
CKC 2 / 3 k~~ 5 / 3 , where e is the average energy dissipation 
rate per unit mass and cr is the universal Kolmogorov 
constant. While experiments have found the effects of in- 
termittency to be important for high-order moments, the 
correction to the spectral form is quite small. Further- 
more, in both experiment [3] and direct numerical sim- 
ulations [H [5] the velocity in homogenous and isotropic 
turbulence is found to exhibit near-Gaussian statistics. 

Quantum fluids, however, are typically described as 
a mixture of two interpenetrating fluids [BJ, a viscous 
normal fluid and an inviscid superfluid exhibiting long- 
range quantum order. There is no conventional viscous 
dissipation in the superfluid component and vorticity 
is confined to one-dimensional quantized vortices which 
possess circulation values that are integer multiples of 
k = h/m — 9.97 x 1CP 4 cm 2 /s [6j. Thus, quantum tur- 
bulence takes the form of a complex tangle of atomically- 
thin vortex filaments of quantized strength [7] . Dissipa- 
tion in the superfluid component for 1.70 K < T < 2.05 
K is mainly produced by mutual friction [5] between the 
quantized vortices and normal fluid. 

Despite these fundamental differences, there have been 
notable studies demonstrating similarities between quan- 
tum and classical turbulence [HI CTJ [TH HH [THl EH [Til [M 
IT?! ITS] . Even though the quality of the supporting ev- 
idence has been questioned 19J, it will be summarized 
here. Experiments on turbulence generated in 4 He by 
two counter-rotating disks observed Kolmogorov energy 
spectra that were indistinguishable above and below the 
superfluid transition [TJ]. The Kolmogorov energy spec- 
trum was seen in numerical simulations of the Gross- 
Pitaevskii equation with small-scale dissipation added to 
the otherwise energy-conserving dynamics |15j . The clas- 
sical decay of vorticity [TO] has been observed in towed 




FIG. 1: Depiction of one-dimensional topological vortices that 
reconnect by merging at the moment to and then separate [25] ; 
the minimum separation distance is S(t) with 8(to) = 0. 



grid PJjl [Tj)] , thermal counterflow [T7] , and impulsive spin 
down [18 experiments. In all these studies the flow scales 
observed were considerably larger than typical intervor- 
tex spacings. These results may be attributed to the fact 
that on such scales the pairwise interactions of quantized 
vortices arc insignificant while the normal and superfluid 
components become "locked" as a result of mutual fric- 
tion. 

In this Letter, we study the velocity statistics of 
quantum turbulence generated by a thermal counter- 
flow on length scales between our experimental resolu- 
tion (^1 fim) and the typical intervortex spacing (~0.1-1 
mm) [17]. On such length scales the interactions of in- 
dividual quantized vortices are important. Specifically, 
quantized vortex reconnection [5DJ [5TJ [55] [23] [H] , where 
two vortices merge at a point, change topology by ex- 
changing parts, and separate (Fig. 1 and [25]), produces 
high, atypical velocities. By analyzing the trajectories 
of micron-sized solid hydrogen tracers, we may compute 
both the velocity statistics of quantum turbulence, and 
identify and assess the effects of individual reconnection 
events. Previous studies have shown that hydrogen trac- 
ers can be trapped on quantized vortices or, if not near 
a vortex, move with the normal fluid under the influence 
of Stokes drag [26] [27] [28] . 

Our experiments are conducted in a cylindrical cryo- 
stat of 4.5 cm diameter using liquid 4 He. The long axis 
of the channel is vertical with four 1.5 cm windows sepa- 
rated by 90°. Particles are produced by injecting a mix- 
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FIG. 2: Time- varying pulsed counterflow velocity distribu- 
tions at T = 1.90 K showing (a) v x and (b) v z for a portion 
of a heat pulse of 0.17 W/cm 2 with the heater turned off at 
t — toff [30|. White denotes amplitudes with zero probability. 



ture of 2% H 2 and 98% 4 He into the liquid helium above 
the superfluid transition temperature (2.17 K) [57]. The 
volume fraction of hydrogen is ~ 10~ 8 — 10 _7 , which re- 
sults in each evident vortex having only a few trapped 
particles so that their effects may be neglected [28] . The 
fluid is then evaporatively cooled to the desired temper- 
ature in the range 1.70 K < T < 2.05 K. The hydrogen 
particles are illuminated by an argon-ion-laser sheet that 
is 8 mm tall and 100 /im wide. Optical laser power varies 
between 2 and 6 W. A video camera gathers 90° scattered 
light with a resolution of 16 /urn per pixel at 60, 80, or 
100 frames per second. 

We study the velocity statistics in decaying quantum 
turbulence initiated by a reproducible thermal counter- 
flow. A spiral nichrome wire heater at the bottom of 
the channel 7.5 cm below the observation volume drives 
the counterflow. The fixed heat flux varying from 0.064 
to 0.17 W/cm 2 drives the system for approximately 5 s, 
after which it relaxes for about 10 s before the process 
is repeated. A 2D particle-tracking algorithm with sub- 
pixel precision [25] extracts single particle trajectories. 

We characterize the resulting dynamics by analyzing 
the particle trajectories |28j . Time- varying distributions 
of the horizontal and vertical velocity components, v x 
and v z , computed by forward differences, are shown in 
Fig. 2 for a typical thermal pulse [515]. The v x distribu- 



tions are always peaked near zero. However, as predicted 
by the two-fluid model [6], the v z distributions exhibit a 
different behavior, since entropy injected by the heater 
is carried upward (v z > 0) by the motion of the nor- 
mal fluid. To conserve mass, the superfluid component 
moves downward opposing the normal fluid motion. The 
bimodal v z distributions when the heater is on represent 
particles with v z > moving upward with the normal 
fluid owing to Stokes drag while particles with v z < are 
trapped in the vortex tangle that moves downward. Once 
the heat pulse ends, the vertical velocities collapse to dis- 
tributions peaked near zero. The probability distribution 
functions (pdf) of velocity components derived from all 
trajectories for times after the heat pulse is turned off 
(i > t s) in Fig. 2 are shown in Fig. 3(a). We focus on the 
tails of these distributions, which are composed of trajec- 
tories with high, atypical velocities; these we attribute to 
quantized vortex reconnection as explained below. 

Near the reconnection moment, quantized vortices 
move with velocities much higher than the background 
flow [25]. These reconnections have been experimentally 
visualized using hydrogen particles |24j and studied nu- 
merically [2U1 |2TJ [22] and analytically [25] . As discussed 
in some of these works, the minimum separation dis- 
tance between reconnecting quantized vortices S(t) (Fig. 
1) evolves approximately as a square-root in time, i.e. 
S(t) = A-J n\t — trj|> where to is the reconnection moment 
and A is a dimcnsionlcss factor of order unity. Thus, for 
lengths between the vortex core radius and the typical 
intcrvortex spacing we expect the velocity to scale as 



v(t) cx \t~t \-^ 2 , 



(1) 



which grows much larger than typical fluid velocities 
when t — *■ t (although cut off by the speed of sound). 

In the pulsed counterflow experiments, a reconnection 
event is evidenced by a pair of nearby tracers rapidly 
approaching or separating. Given the large number of 
possible particle pairs analyzed (~ 10 10 ) an ad hoc cri- 
terion is needed to select likely reconnection events. We 
define a pair of particles i and j as marking reconnection 
at time t if the separation 5ij(t) = |r ^ (i ) — Tj(t)\ satisfies 



S ij (t± 0.25 s)/S ij (t) >4, 



(2) 



where Vi(t) is the two-dimensional projection of the po- 
sition of particle i at time t and the plus (minus) sign de- 
notes particles that separated after (approached before) 
an event, which we label as forward (reversed) events. 
The duration of 0.25 s is chosen since greater times are 
dominated by boundary effects and the presence of other 
vortices. The criterion ([2| excludes all but a fraction of 
possible pairs leaving ~4x 10 4 reconnection events. 

The measured separations 5(t) for four typical forward 
events are shown in Fig. 4(a). Solid line fits invoking a 
correction factor to the predicted scaling as 



S Rt (t) = A[ K (t-t )} 1/2 [l + c(t-to)l 



(3) 
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FIG. 3: Local velocity and energy statistics derived from the 
data in Fig. 2 for all particle trajectories with t > t a g (com- 
puted from over 1.1 x 10 6 values of velocity). All distributions 
are scaled to give unit variance using a Vj . = 0.066, o v , = 0.074 
cm/s or ob = 0.017 (cm/s) 2 . (a) Probability distribution 
function of v x (black circles) and v z (red squares). The solid 
lines are fits to Pr„(«i) = a\vi — Wi~ 3 , where i is either x 
(black) or z (red) and vl is the mean value of Vi. For compar- 
ison the dashed (blue) line shows the distribution for classical 
turbulence in water [32] computed from over 10 7 velocity val- 
ues. The distribution is scaled using a v — 0.25 cm/s, and with 
a peak value matched to the v x data. The velocity statistics 
in water are close-to-Gaussian over five decades in probabil- 
ity, (b) Probability distribution function of E = (v x + t> 2 )/2 
with a fit of the form Pte(E) = aE~ 2 shown (in red). 



describe the data well [31 . The fits minimize \ 2 = 
n ~ 1 S™=i [(^f t— ) / cr] 2 , where i denotes the movie frame, 
a — 4 /im (0.25 pixels) is an estimate of the uncertainty 
of the particle positions, and n = 15, 20, 25 for data col- 
lected at 60, 80, or 100 frames per second, respectively. 
To fit reversed events, we use the same form ^ with 
(t — 1 ) replaced by (t — t). The distributions of A and c 
for forward (reversed) events, determined from fifty heat 
pulses, are shown in black (red) in Figs. 4(b) and (c). 
These are calculated only from events with \ 2 < 4; about 
50% of the pairs satisfying (|2J) meet this \ 2 criterion. 
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FIG. 4: Statistics of the reconnection fits, (a) Four typical 
forward events. Symbols denote the measured separation 8(t) 
of pairs of particles on reconnecting vortices with an example 
error bar a — 4 fim while solid lines show fits to ^3J. (b) 
Normalized distributions of the amplitude A for 20,300 for- 
ward events (black circles) and 19,600 reversed events (red 
squares), (c) Normalized distributions of the correction am- 
plitude c for the same 20,300 forward events (black circles) 
and 19,600 reversed events (red squares). 



To model the pdf of the velocity derived from particle 
trajectories, we may use the transformation 



Pr„(«) = Pr t [t(v)]\dt/dv\ 



(4) 



where Pr v (v)dv is the probability of observing a velocity 
between v and v + dv at any time while Pr t (t) dt is the 
uniform probability of taking a measurement at a time 
between t and t+dt. Hence, accepting the scaling relation 
([lj, we predict for large v (small t) the behavior 



Pr^(w) oc \dt/dv \ oc \ v 



-3 



(5) 



The v x and v z pdfs derived from all particle trajecto- 
ries for t > t Q g for the same pulse in Fig. 2 are shown 
in Fig. 3(a). The solid lines are fits to ^ allowing for 
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a mean flow. To emphasize the distinction with clas- 
sical turbulence, a velocity pdf from an oscillating-grid 
experiment in water [32] is also shown. Evidently the ve- 
locity pdfs in superfluid helium differ drastically from the 
near-Gaussian velocity pdfs observed experimentally [3], 
and in direct numerical simulations of homogenous and 
isotropic classical turbulence [H \5\ . One must note, how- 
ever, that tracer particles in superfluid helium respond 
separately to the normal fluid and the quantized vortices 
(which are influenced by the normal fluid and superfluid) ; 
this is fundamentally different from that in water. 

By the same argument, the tails of the pdf for the 
kinetic energy per unit mass E — (v% + v1)/2 will be 
dominated by reconnections. Accepting the relation (JT|), 
we have E(t) oc \t — | 1 and so for large E we expect 



Yx E {E) oc \ dt/dE\ cx E~ 



(0) 



The pdf of E computed from the data in Fig. 3(a) (which 
includes all particle trajectories) is shown in Fig. 3(b). 
The departure from the predicted power-law behavior for 
low velocities and energies may reasonably be attributed 
to effects from the boundaries and nearby vortices as well 
as to the background drift of the normal fluid. 

In conclusion we have shown that the velocity statistics 
of quantum turbulence in superfluid 4 He differ drastically 
from those for classical turbulence owing to the topologi- 
cal interactions of vortices that are different from those in 
classical fluids. Previous studies argued that the interac- 
tions of magnetic field lines can cause the velocity statis- 
tics of magnetohydrodynamic (MHD) turbulence also to 
differ from classical turbulence. The power-law tails in 
the distributions of electron energies observed in astro- 
physical settings (Fig. 3 in [33] and Fig. 2 in [34]) have 
been attributed to magnetic reconnection |35j . Further- 
more, theories for MHD turbulence propose that frac- 
tional diffusion may be the dominant transport mecha- 
nism 3G . Such diffusion is associated with power-law 
tails in velocity distribution functions. Since reconnec- 
tion is a principal dissipative mechanism in supcrfluids 
near absolute zero, superfluid experiments might provide 
another "laboratory" for studying strong MHD turbu- 
lence, as well as other systems exhibiting one-dimensional 
topological defects [37] such as liquid crystals, supercon- 
ductors, Bose-Einstein condensates, and cosmic strings. 
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